An alternative construction of the positive inner product for pseudo-Hermitian 
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This paper builds on our earlier proposal for construction of a positive inner product for pseudo- 
Hermitian Hamiltonians and we give several examples to clarify our method. We show through 
the example of the harmonic oscillator how our construction applies equally well to Hermitian 
Hamiltonians which form a subset of pseudo-Hermitian systems. For finite dimensional pseudo- 
Hermitian matrix Hamiltonians we construct the positive inner product (in the case of 2 x 2 matrices 
for both real as well as complex eigenvalues). When the quantum mechanical system cannot be 
diagonalized exactly, our construction can be carried out perturbatively and we develop the general 
formalism for such a perturbative calculation systematically (for real eigenvalues). We illustrate 
how this general formalism works out in practice by calculating the inner product for a couple of 
VT symmetric quantum mechanical theories. 

PACS numbers: 03.65.-w, 03.65.Sq 



I. INTRODUCTION 

In a recent paper ^ we described a systematic pro- 
cedure for constructing the positive inner product for 
a quantum mechanical systems described by a pseudo- 
Hermitian Hamiltonian which satisfies [2] 



h = s-^h''s. 



(1) 



Here 5 is a bounded operator which can be chosen to 
be Hermitian (in the sense of Dirac). The crucial con- 
cept in our construction is the knowledge of the genera- 
tors of energy eigenstates which acting on a given refer- 
ence state generate all the eigenstates of the Hamiltonian. 
The method works equally well for systems with real or 
complex energy eigenvalues. We note that a pseudo- 
Hermitian Hamiltonian reduces to a Hermitian Hamil- 
tonian when S — 1 and, therefore, our construction of 
the positive inner product also reduces to the standard 
Dirac inner product in this case (which we discuss in 
section III). In [l| we illustrated the method through 
the example of the Lee model [3, |3| with an imaginary 
coupling constant [sj. In this paper we give additional 
examples within the context of finite dimensional ma- 
trix Hamiltonians (which are pseudo-Hermitian) to clar- 
ify our method. 

However, most quantum mechanical systems cannot be 
solved exactly. As a result, it is not possible to determine 
the energy eigenstates of the system and, therefore jtheir 
generators in a closed form. As we observed in [l|, in 
such a case, the generator as well as the inner product 
can be determined only perturbatively. (The pseudo- 
Hermitian systems differ from the Hermitian systems in 
this respect, namely, the inner product of the system 
depends on the dynamics of the system in a nontrivial 



manner 

0, i, S iO In this paper, we would like to 
develop the idea of a perturbative determination of the 
inner product in detail for real energy eigenvalues (the 
discussion for complex energy values is straightforward) 
and work out various VT symmetric 0] examples to 
illustrate the method. 

The paper is organized as follows. In section II, we 
briefly recapitulate our proposal for constructing the pos- 
itive inner product in a pseudo-Hermitian quantum me- 
chanical system, clarifying some of the details not fully 
explained in In section III, we apply our method to 
the case of the harmonic oscillator which can be taken as 
the unperturbed Hamiltonian in some of the examples we 
discuss. In section IV, we work out in detail a pseudo- 
Hermitian 2x2 matrix model with real and complex 
eigenvalues in an analogy with Q. We also discuss the 
method in the context of n x n matrix models. In section 
V we describe the general formalism for the perturba- 
tive determination of the generator of energy eigenstates 
and, therefore, the perturbative construction of the posi- 
tive inner product (for real energy values). In section VI, 
we apply the general method to various VT symmetric 
quantum mechanical models with real energy eigenval- 
ues. We conclude with a brief summary in section VII. 



II. RECAPITULATION OF THE 
CONSTRUCTION OF THE INNER PRODUCT 



In this section, we briefly recapitulate the essential 
points discussed in [1] in the construction of the inner 
product in a pseudo-Hermitian quantum mechanical sys- 
tem described by ([T|). We showed that if we define an 
operator 
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q = SA, 



(2) 
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where [A, H] = 0, and choose A properly, then the 
quadratic form 



(3) 



defines an inner product (on a suitably defined Hilbert 
space, see [l| for details) satisfying 



(4) 



with a unitary time evolution. Here and in what follows 
a bar denotes complex conjugation. The operator q can 
be systematically constructed from a knowledge of the 
generators of the energy eigenstates of the theory in the 
following way. 

An operator (7^ satisfying the relation 



H(Te = E(7e + CTb/Ce, 



(5) 



is defined to be a generator of the eigenstates of H with 
eigenvalue E provided 

{i) there exists at least one vector |-!/)), solving 

fcB|?/')=0, for all £; e spect(i/), (6) 
with cr^; !■(/;) ^ 0. 
(m) there exists at least one vector |0) solving 

4|0>=O, for all e spect(i/), (7) 

with (i/jI^) ^ 0. 

{Hi) Ue has an inverse Oe^, at least acting on the refer- 
ence state l-i/)), such that 



CTe VeIt/)) = IV'), 



(8) 



and an adjoint inverse, (ctb) ^, well-defined when 
acting on |(/)), with \4>) ^ 0. 



From (O and ©, it follows that 

IV'e) = CTeIV'), 



(9) 



is an eigenvector of H with eigenvalue E. Similarly, from 
the (Dirac) adjoint of ([5]) as well as using ([7]), we can 
show that 



(10) 



is an eigenvector of with eigenvalue E. Furthermore, 
if Pe denotes projection operator on to the energy eigen- 
state IV'e) with energy E 



-PeIV'e') = &EE' IV'e), 



(11) 



then we showed in ^ that the action of q can be ex- 
pressed as 



q{aEm = (4) 



t^-l| 



As a result q has the operator form 

9 = E('^^)"''7ot-1Pe, (13) 

E 

where we have identified 

10) = golV), such that = (V'kolV') = 1- (14) 

We note that in practical calculations, it is sufficient to 
have a partial inverse [*1 cr"-'^ as we will demonstrate in 
the next section. 

The construction ([5|)- P^ is quite general. It holds 
for pseudo-Hermitian Hamiltonians with real or complex 
eigenvalues. It also works for Hermitian Hamiltonians 
(which is a subset of pseudo-Hermitian Hamiltonians) for 
which S = \ and in this case the operator q can be chosen 
to be g = 1 leading us back to the Dirac inner product. 
In [Ij] we applied this construction to calculate q for the 
Lee model [3, 01 with an imaginary coupling The 
purpose of the present paper is to work out more exam- 
ples of diverse nature in order to shed more light on this 
method. It was pointed out in [l| that when the exact 
eigenstates of the Hamiltonian are not known, the gener- 
ators as well as q can be constructed perturbatively and 
one of the goals of the present work is to describe sys- 
tematically how such a perturbative calculation works. 
In the next section we apply our method to the case of 
the harmonic oscillator which can be taken as the zeroth 
order Hamiltonian in some of the subsequent examples 
where we determine q perturbatively. The Hamiltonian 
for the harmonic oscillator is, of course, Hermitian and 
this would also show how our method leads back to the 
standard Dirac inner product in cases where the Hamil- 
tonian is Hermitian. 



III. HARMONIC OSCILLATOR 

Before we determine q for finite dimensional matrix 
models and discuss how q can be determined perturba- 
tively in systems where the exact energy eigenstates are 
not known, let us describe how the construction of the 
last section works in the simple example of the harmonic 
oscillator. In this case, we expect the inner product to co- 
incide with the standard Dirac inner product and, there- 
fore, we expect to be able to show that we can choose 

The Hamiltonian for the harmonic oscillator can be 
written as 



H\i.o. — 



1) 



a) a 



Hi 



(15) 



(12) 



[*]A valid partial inverse can be found by noting that the generators 
are not unique. In fact, if cte is a generator for H then cte = 
f{H)crE is also a generator (with the same reference state I?/")). For 
the harmonic oscillator we get that (tT~^)t = 1/a as the generator 
for , but we prefer to use instead //^((j— l)t = a^. 
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where a, denote respectively the annihilation and the 
creation operators defined by 

a= ^{x + ip), ^ -^{x~ip), (16) 



and we have subtracted out the zero point energy in ([T| 
for simplicity. Here we have also set m = = 1 = fi, for 
simplicity. The energy eigenstates of the oscillator are 
well known and satisfy 

Hh.o.\4'Ej ^E^\i;^J ^n\^^), n = 0,l,2,---. (17) 

Therefore, we can write the normalized eigenstates as 



(at) 



l^o), a|V'o)=0, (18) 



where \iIjq) denotes the ground state of the system (com- 
monly denoted by |0)) and the energy eigenstates satisfy 
the orthonormality relation (with respect to the Dirac 
inner product) 



71/ — ^mn- 



(19) 



We see from p8)) that we can identify the reference 
state with the ground state and the generator of the en- 
ergy eigenstates as 



IV-) = IV'o) 



CTr„ — (Jr, 



(at) 



(20) 



Furthermore, from the canonical commutation relation 
[a, at] = 1, it follows that 



Hh.o.CTn — EnCTn + (JnHh.o., 



(21) 



where we have used En = n. Comparing this with ([5]) wc 
determine that in the present case 

^ = -ffh.o.j (22) 

which indeed satisfies 

fcnIV') =i?h.o.|^o) =0. (23) 
It also follows that for the present problem 

\^)^\^P) = \^o), 90 = 1. (24) 

The partial inverse of the generator ([8]) can also be 
determined from the canonical commutation relation to 
be 



(25) 



We emphasize here that only a partial inverse is necessary 
for our construction which allows us to avoid awkward 
terms of the form (at)"" for the inverse. We also note 



that the projection operators for the energy eigenstates, 
in this case, can be written in the simple form 



n 

which satisfies 



{a 



P ^ P 



(26) 



(27) 



tn\ \fn\ 

With all these, we can now determine q (see (fT3|)) to be 

n 

SO that the inner product ([3]), in this case, coincides with 
the standard Dirac inner product 



(29) 



It is worth emphasizing here that the relation ([5]) defin- 
ing the generator does not determine its scale uniquely, 
namely, the generator is defined only up to a multiplica- 
tive factor. If we had defined the generators in ([^0]) as 
On = c„(at)" instead (giving = ^^), (HSl) would 
lead to 



(30) 



Thus, we see that we can choose any coefficient multi- 
plying the generators and it will change g by a constant 
factor at each P„. However, any choice of (t„ enforces 
(V'nlV'n)? = 1. In particular, this determines that g = 1 
for c„ = reducing the inner product to the Dirac 
product. 



IV. PSEUDO-HERMITIAN MATRIX 
HAMILTONIANS 

In this section we discuss in detail two examples where 
we describe how our method applies to finite dimen- 
sional matrix Hamiltonians which are in general pseudo- 
Hermitian. The finite dimensional matrix Hamiltonians 
can be solved exactly in principle. 



A. 2x2 matrix 

As the first example, we apply our method to a sim- 
ple 2x2 matrix Hamiltonian with real as well as com- 
plex eigenvalues. We recall that the T^T symmetric 2x2 
matrix model described by the Hamiltonian (see for 
details) 



re** s 
s re' 



(31) 
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has been studied extensively in the past. Let us, there- 
fore, analyze the generalized 2x2 matrix Hamiltonain 
(r, s, t, 6, (j) are real parameters) 



H 



r e 
te-'' 



s e 
r e 



(32) 



which is not Herniitian unless s = t,6 = Q. For s = t,(j) = 
this model reduces to the VT symmetric theory (|3T 
where the parity operation is identified with 



1 

1 



(33) 



with T denoting complex conjugation. The Hamiltonian 
in (|32p. however, is not VT symmetric in this context. 
On the other hand, if we define a "generalized parity" 
operation through the 2x2 matrix 



V = 




(34) 



and choose time reversal to correspond to complex con- 
jugation, then it is easily verified that the general Hamil- 
tonian (|32|) is VT symmetric, namely. 



VTHVT = H. 



(35) 



The "generalized parity" operator in ((34| reduces to 
when s = t. In general, however, let us note that 
while this operator is idempotent, it is not Hermitian, 



1, v'f^v. 



(36) 



Furthermore, unlike the case in VT symmetric theories 
d, 01 (where VHV = and, therefore V defines the 
matrix 5' in ([l])), here we have 



VHV = H* ^ ijt. 



(37) 



Therefore, V cannot be identified with S* in ([T]) which 
is used in our construction. However, we can define an 
operator (2x2 matrix) 



S 



which leads to 







J4> 







(38) 



(39) 



In other words, even though P does not take the Hamil- 
tonian to its Hermitian conjugate, S does and the gen- 
eral Hamiltonian H in Ij32p is pseudo-Hermitian and our 
method can be directly applied. 

The energy eigenvalues of this system are given by 



E± 



i9 ± Vst—r'^ sin^ ( 



(40) 



and they are real for st > sin^ 6 while they are complex 
for st < sin^ 9. (We do not consider the degenerate 
case st = sin^ 0, for which the Hamiltonian cannot be 
diagonalized, simply because it would introduce pseudo- 
eigenvectors.) We would analyze the two cases separately 
in the following. 



1. Real eigenvalues: 



In the case of real eigenvalues, let us define 



Q = \J st - r2 sin^ B = real, (41) 

so that the two real energy eigenvalues can be written as 

E±=r co'&B±Q = E±. (42) 

The corresponding eigenvalues can be determined to have 
the forms 

_ ^_ ( (if' VQ + ^rsmOe^^/^ \ 
1^^+^^ /— r-rl /,^l/4 ^ a ^-»0/2 ) ' 



ir sm t) e 



ilfy Q-zr sine e ^^/' ] 

which have been normalized in the conventional sense 
for simplicity (although it is not necessary). We note 
that the two energy eigenstates are also eigenstates of 
VT with eigenvalue 1 (they correspond to singlet states 
under VT), namely, 

VT\ib^^) = \^,^), VT\^^_) = \ij,_). (44) 

The projection operators onto the two states are given 
by 

_ J_ fQ + ir sine se"*' 
+ ^ 2Q\ te-"t' Q-ir sin( 

^ J_ (Q^ IT sine ~se^* \ 

2Q\ -te-'^ Q + ir sine)- ^ ' 

Let us choose the reference state (see ([SJ or ([S])) 

IV) = IV-J. (46) 

It follows now from (|43p that the generators of the states 
are given by 

/ 




= 1, k^_= -k^^=2QP^^. (47) 
The reference state ([7|) can also be determined to be 



2Q 



2Q 



(48) 

where the reference state has been normalized according 
to p4| . It follows now from (fT3|) (using (|47p as well as 
\) that 



s + t 



t 



—ir smt/e 
2Q2 yir sinBe-'-'f' s 



(49) 



It can be checked now that 

(V'eJq'IV'e,) = (V'eJV'e,)? = (^ij, j,j = ±- (50) 
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2. Complex eigenvalues: 



where 



For the case st < sin^ 6, as we have noted the eigen- 
values are complex and let us define 



\/ sin^ 9 — st — real, 



(51) 



so that the two energy eigenvalues ((40|) which become 
complex conjugate pairs can be written as 

E = r cos9-iQ, E cos9 + iQ. (52) 

The two eigenstates can be determined to have the 
forms 



(-0 



(s + t)r sm9 + {s - t)Q 




iip/2 



iy s{r sm9 — Q)e 

t{r sm9 + Q)e~"f'^^ 
i 

{s + t)r sin6' + (s - t)Q 
s{r sin6' + Q)e''^/^ 



X I • r , I ' (53) 

'i\Jt{r sin6'- Q)e""^/2 

It s clear that the eigenstates (|53|) of the Hamiltonian for 
complex eigenvalues correspond to a doublet representa- 
tion of VT (compare with (|44p ). namely, 



(54) 



The projection operators onto these two states take the 
forms 

p _ J_ f-{r sine* - Q) ise">' 
2Q\ ite''-'^ r siT\9 + ^ 

_ ^(r^\n9 + Q -zse^* \ 
^ " 2Q I -^te-^^ -{rsm9-Q)J- ^''^ 

Let us next choose the reference state 

W = (56) 

It follows now from (|53p that the generators of states are 
given by 







(Jr = 1, kr: 



k^ — ^^^QPe- 



(57) 



The reference state ^ can now be obtained to have the 
form 

i,\ I ,\ {s + t)r am9 + {s-t)Q . 



{s + t)r sin 6* + {s ~~ t)Q 







-Vh 



(59) 



2^/s<Q 

It can be checked that this state satisfies 

410) = 4|(/.) =0, = (^|go|^) = 1. (60) 

It follows now from (using ([Ff|) as well as ([55]) ) 
that 



1 = -CTeIoPe - QoCTePe 

{s + t)r sm9^ (s - t)Q 



2VstQ 
(s + t)r sa\9+ [s - t)Q 



S{Pe+Pe) 
S. (61) 



2y/stQ 

It is easy to check that 

{^E\^E)g^O^ {^e\^e),, (VfilV^), = 1. (62) 

In this case, we note that the operator q simply scales the 
operator S in ([55)1 which leads to the pseudo-Hermitian 
nature of the Hamiltonian. 

B. n X n matrix 

Here we will solve for the operator q for an n- 
dimensional Hilbert space, so that q will be an rt x n 
matrix. We will assume that the Hamiltonian H has n 
real, distinct eigenvalues (no degeneracies) so that H has 
a diagonal Jordan decomposition. Therefore, we can find 
an invertible matrix R such that 



R ^ HR — Hoi 



(63) 



where Hq = diag{i?i, i?2, En}- Because Hq is diagonal 
there exists a basis of vectors, \vi) satisfying, 

Ha\v,)=E,\v,), {vj\v,)^5y,, E, = E,, (64) 

for 1 < i, j < n. 

However, in the non-diagonal basis, the eigenstates for 
H are given by {ipi) = R\vi). The Dirac-inner product of 
two such states is 



.\R^R\v,) 



(65) 



and is generally not orthonormal with respect to the 
eigenstates We can modify the inner product by 

introducing an operator q as discussed in the previous 
section. 



(66) 



2VstQ 



In order to restore orthonormality it would appear from 
((65| and ^ that we should choose q = {R^^YR^^. On 
the other hand, we recall from (fT3l) that 



(67) 



9 = '^i<^El)^90'TElPi, 
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where ct.eJV') = IV'j) defines the generators of states, 
Pj\4^i) — Sijlipi) defines the projection operators, and 
Qo is the matrix satisfying go IV') = (V'l'/') — for the 
reference states 

To show that the two definitions are in fact equivalent, 
let us choose the reference state to be — \ipn)- Since 
\tpn) — ^l^'n), this imphes that \(f>) — (i?^^)^|w„) (be- 
cause we want = 1 to hold). Furthermore, since 
golV'n) = golV") = 10), then using P„ = R\vn){vn\R^^ = 
\'4'n) {ipn\q we can write the identity 



qoPn = qPn- 



(68) 



(We note that although qo is generally different from q, 
it must act the same way on the reference state |'0n)-) 

Second, we observe that in the diagonal basis, \vi), 
there exists an upper-triangular matrix t satisfying 
t\vi) = \vi+i), with t\vn) = \vi), and — t^^ [t is just 
the matrix with ones above the diagonal, and one in the 
lower left-hand corner). This gives = cr* = Rt^R^^ 
(here denotes the i th power of t), where a — RtR^^ 
defines the generator of states. The statement a^'^Pi = 
Pn<y~^ is evident. 

Using our definitions above, we can calculate directly 

n 

n n 
= Y.i^~yq^PnO~' ^Y.^o-yqP^o-^ 

n 

= J2^{R~'yt'R^){{R-yR-^)Pr,{Rt-''R-^) 

i=l 



= {R-yR-\ 



(69) 



In finite dimensional Hilbert spaces, such as this, it is 
probably easier to calculate q = {R~^y R'^ directly from 
the Jordan matrices R (which consist of each eigenvector 
of H in one of the columns). However, in infinite dimen- 
sional Hilbert spaces there is no Jordan decomposition, 
but the operators ctb which generate the eigenvectors do 
exist. Thus, our expression for q is valid in both finite, 
and infinite dimensions, whether or not one can construct 
the Jordan matrix. As we will see in the next two sec- 
tions, there exist operators (also denoted by R) which 
will play a role very similar to the Jordan matrix given 
above. These will turn out to be the perturbation oper- 
ators. 



V. PERTURBATIVE DETERMINATION OF q 

Let us next consider a pseudo-Hermitian system with 
real energy eigenvalues (as in VT symmetric theories 
[1, 0) where the energy eigenstates may be difficult to 



determine exactly. In this case, let us assume that we 
can write the pseudo-Hermitian Hamiltonian H as 



H = Hn 



eV{x), 



(70) 



where Hq is the part of the Hamiltonian H which we can 
diagonalize and we treat eV{x) as a perturbation. Of 
course, Hq can be Hermitian or non-Hermitian. How- 
ever, in most practical examples, it can be chosen to be 
Hermitian and this is the case we will discuss here. Let 
us denote the generators for the eigenstates of Hq by CTe"^ 
so that we can express the eigenstates of Hq as 



1^ 



(0)\ 



(71) 



where denotes the reference state for the system (see 
(|n]) and Q). Since Hq is chosen to be Hermitian, as in 
the case of the harmonic oscillator, it follows that (see 

(El) 



), qo = 1. 



(72) 



Thus, for the (diagonalizable) unperturbed part of the 
Hamiltonian we can determine (see (f28| ) 



-,(0) 



5:(af -^)t,o(4"^)-^pi°^ 

E 

J2{aT -y ia^^Y' Pi"' = 1, (73) 



as we have seen in ([28|). Here P^"^ = I'Vi"^) (V'e''' I is the 
projection onto an eigenstate of Hq with eigenvalue 

The eigenstates of the total Hamiltonian H in ([70]) 
can be written as a series in powers of the perturbation 
parameter e, namely, (e* denotes the ith power of e) 



(74) 



where IV'e'') denotes the ith order correction to the state 
IV'^^). We can always choose the corrections to a given 
state to be orthogonal to the state, namely. 



(Vi"Vi*^) = o, *>o. 



(75) 



In this case, the perturbative corrections to the state 
IV'i*^^) can be obtained to have the compact recursive form 
(i > and we are assuming the energy eigenvalues to be 
discrete) 



1^ 



i=i 



(76) 
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where i?*^*^ denotes the operator that transforms the un- 
perturbed state IV'is'^) to IV'e'') and 



(77) 



From (|76p we note that we can define the correction to 
the generator of the state at the ith order as 



(78) 



with the same (zeroth order) reference state \tp). There- 
fore, to any order in perturbation we can write 



R 



1=0 



(79) 



i=0 



where i?'-*^-' = 1. By taking the adjoint of ([70)1 . and using 
|(/)) = (see (l72l) '). we likewise determine the corrections 
to the eigenstates of H^^ to any order to be (see PH)) ) 

10,) - (i?-i)t|(/)i")> = {R-y{<7^°^ IV), (80) 

with the same reference state | ■(/'). 

Puting this into the equation for q (see (|13p ) we find 
that for real (recall that we are considering a pseudo- 
Hermitian system with real energy eigenvalues as in VT 
symmetric theories) 

E 

= (i?-i)t 9(0) (81) 

where q^^^ is the q operator for the unperturbed Hamil- 
tonian which, for a Hermitian ifo, is given by q^'^^ — 1 
(see ([75)) '). In this case, therefore, ([5T|) reduces to 



(82) 



For real eigenvalues, we see that q only depends on the 
perturbation operator i?, and that the problem is solved 
once we determine R from perturbation theory. If we had 
not chosen iJo to be Hermitian, then g^^^ in ([STjl would 
be nontrivial which would need to be determined as well. 



VI. EXAMPLES 

In this section we would work out two examples of VT 
symmetric theories (with real energy eigenvalues) to il- 
lustrate how the perturbative calculation is carried out in 
practice. We recall from ((82|) that when the unperturbed 
Hamiltonian is Hermitian, then to any order in the 
perturbing parameter, the operator q is determined from 
the operator R defined in ([75)1 - (1751) to that order. 



A. Example: H ^ \{p^ + - I) + iex^ 
Let us consider the Hamiltonian given by p^ . 

H = i(p2 iex^ = Ho + eV{x), 



(83) 



where e is a real constant parameter and we recognize 
that Hq — i?h.o. describes the harmonic oscillator which 
we have studied in section III (we continue to identify 
m ^ Lu = I = h). This model is known to be VT sym- 
metric with real energy eigenvalues and has been studied 
extensively in the past. Here we would calculate the op- 
erator q (which defines the inner product) associated with 
this system perturbatively using our method. 

The Hamiltonian for the harmonic oscillator is Hermi- 
tian and the system can be solved exactly. As a result, 
we can choose Hq = i/h.o. to be the unperturbed Hamil- 
tonian in which case we can take over the analysis of 
section V. We have already seen in (fT5)l and (|20p that 
the generator for the harmonic oscillator has a very sim- 
ple form. Therefore, for the unperturbed system we can 
identify 

l^i"^>=4°^IV^o)=af IV'), i?(°)=n. (84) 

Furthermore, recalling that x — -^(a + a^) (see, (fT6|)). 
the first order correction to the unperturbed eigenstate 
can be determined from (|76p as 



- p(0) _ (0) 1^™/ 

= ^{~1{a~a^f + {ia + a^na-a^)}) 



|p3 + a;2p-zx) |Vi°)) 



-p" + x~p — tx \ cr^°) I tp) ■ 



(85) 



Thus, comparing with (|79p we determine that to first 
order in perturbation 



R = i?^"^ + ei?^^^ = 1 - e ( ^p^ + x^p - ix 



(86) 



Similarly, carrying out the perturbative calculation to 
order we determine (see ([79| ) 



R^l + ei?(^) + e^i?'^' -I- e^i?(^\ 



(87) 
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where 










- x^p 




/23 6 
V288^ 


23 




^96' 



X p 



= ( 



59 



„9 



P 



m 

48 
59 

288 



11 . 



'^4 2 

96 ^ 
41 



2 7 

a; p 



-a; — 
593i 



6 



V1296 

109 4 5 1591i 

288 ^ ~ 288 

307 , 649i n 
x > + — — xp 



288 



288 
769 



163 
120 



x^p' 



864 
17 



48 



48 

443i , 2 



41 



685i 



287i 



X 



72 288 
To this order, therefore, we can calculate 

^3 



-e3(i?(3)_^(^(l)y 



,(89) 



And from this we can calculate q to order (see (|82|) ) 
to have the manifestly Hermitian form. 



g = 1 - e 
137 

I ^3 



4 

px p 



2xpx 

41 177 

a; a 

144 8 

181 41 4 4 

p a; px 

3 ^ 36 ^ 



169 
144^ 



169 

Is 



■xp x 



523 
108 
253 
"36" 



x^p^x^ + 



x^p^x^ 



3463 



5 



-P' . 

60 ^ 
155 7 
36 ^ 



457 
5669 



X px 



xp^x 



72 

^p^)+0(.^). (90) 



This can be compared with ^lOj to find linear order agree- 
ment. The difference in higher order terms is due to our 
choice of normalization {^p\^p)q = — 1- 

It is worth noting here that in this case the perturbing 
Hamiltonian is parity odd while the free harmonic oscil- 
lator Hamiltonian is invariant under parity (each term is 
individually VT invariant). As a result, the zeroth order 
energy eigenstates \ipn'') are parity eigenstates and at ev- 
ery even order \ipn^) would have the same parity as lipn"^) 
while at every odd order \ipn'^) would have the opposite 
parity. This leads to the fact that at odd order (where 
the correction to the energy can become imaginary) 



0, 



2m -I- 1. 



(91) 



The even order corrections to energy, on the other hand, 
are all real. Therefore, 



En 



1=0 



(92) 



is real at any order of perturbation. This is an alternative 
way of understanding the reality of energy in this model. 



B. Example: H = l{p^ + x'^ - 1) + iax + iex^ 

From the previous results it is a simple matter to calcu- 
late q in the case of the VT symmetric theory described 
by 



H = i(p^ -I- x^ - 1) -I- iax + iex^, 



(93) 



where a, e are real constant parameters. First we note 
that with the canonical commutation relations we can 
write (see ((IS])) 

g-ap/2 jj^^ ^ gap/2 ^ jj^ ^ ^ -^^ _ ^2^^_ (g^) 

Therefore, in this case, we can write (|93| 

H = e-"f /2 ff^ ^ gap/2 + c^2/4 _^ .g^3 ^ + j£^3^ (95) 

where we have identified 

Ho - e-"P/2 jj^^ ^ gap/2 ^ g^2/4^ (gg) 

We can take this Hq to be the free Hamiltonian in this 
case for which 

^(0) ^ ^(h.o.) ^ g^2/4^ 1^(0)) ^ (97) 

We can now carry over the perturbative analysis of the 
previous example on this state and determine q (and, 
therefore, the inner product) to any order from 

q = ((ii'e-"P/2)-i)t(^e-«P/2)-i = (^-i)t e"P i?-i,(98) 

with R calculated in the previous example. 

These two examples illustrate how the perturbative de- 
termination works when the exact energy eigenstates are 
difficult to obtain. 



VII. SUMMARY 

In this paper we have tried to build on our earlier pro- 
posal [l| for determining the positive inner product in 
the case of a pseudo-Hermitian Hamiltonian. The crucial 
concept in our method is the generator of energy eigen- 
states. We have clarified various aspects of our proposal 
which were not possible to explain in the earlier paper. 
We have shown through the example of the harmonic os- 
cillator how our proposal reduces to the standard Dirac 
inner product when the Hamiltonian is Hermitian. We 
have given (additional) examples of finite dimensional 
pseudo-Hermitian matrix Hamiltonians to explain how 
our method works in practice. When the quantum me- 
chanical Hamiltonian cannot be diagonalized exactly, the 
generator as well as the inner product need to be con- 
structed perturbatively. We have developed the general 
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formalism for this systematically for real energy eigen- 
values. Furthermore, we have applied the formalism to 
two VT symmetric Hamiltonians to illustrate how it is 
carried out in practice. 
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